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Abstract 



We use QCD sum rules to calculate the hadronic matrix elements govern- 
ing the rare decays B — > K£ + £~ and B — > K*£ + l~ induced by the flavour 
changing neutral current b — > s transition. We also study relations among 
semileptonic and rare B — » K^-*' decay form factors. The analysis of the 
invariant mass distribution of the lepton pair in B -» K(*k+£- and of the 
angular asymmetry in B — > K*£ + £~ provides us with interesting tests of the 
Standard Model and its extensions. 



PACS:11.50.Li, 13.25 
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I. INTRODUCTION 



Rare 5-meson decays induced by the flavour changing neutral current b — > s transition 
represent important channels for testing the Standard Model (SM) and for searching the 
effects of possible new interactions J]J. As a matter of fact, these processes, that in SM 
do not occur in the Born approximation, are particularly sensitive to perturbative QCD 
corrections and to possible higher mass scales and interactions predicted in supersymmetric 
theories, two Higgs doublet and topcolor models, left-right models, etc. Such interactions 
determine the operators and their Wilson coefficients appearing in the low energy AB = 1 
effective Hamiltonian Hw that governs the b —>■ s transition. 

From the experimental point of view, the radiative b — > S7 decay has been observed and 
measured by CLEO II Collaboration both in the inclusive B — > X s i and exclusive B — > K*^y 
modes; the experimental results 

B(b -> sj) = (2.32 ±0.57 ±0.35) 1(T 4 [2] (1.1) 

and 

B(B° -> iT°7) = (4.0 ± 1.7 ± 0.8) 10~ 5 
B(B~ -> K*~i) = (5.7 ± 3.1 ± 1.1) 10~ 5 



[3] (1.2) 



have prompted a number of analyses aimed at restricting the parameter space of various 
extensions of the Standard Model ||]. Similar analyses have also been proposed for the 
transition b — > s£ + £~, that has not been observed, yet []5|; in this case, the invariant dilepton 
mass distribution and the asymmetry of the dilepton angular distribution, together with the 
total decay rate, can be used to study the features of the interaction inducing the decay. 
However, for the exclusive modes such as B — > K£ + £~ and B — > K*£ + £~ one has to face 
the problem of computing the matrix element of Hw between the states B and K, K*, a 
problem related to the nonperturbative sector of QCD. 

For these matrix elements, either specific hadronization models |^J7|] or information from 
two point function QCD sum rules H and from the heavy meson chiral theory M, embedded 
in the vector meson dominance framework, have been used so far. The resulting theoretical 
predictions are characterized by a considerable model dependence; it should be noticed that, 
differently from the case of B — > K* r y 1 where the hadronic matrix element must be computed 
only at one kinematical point, in correspondence to the on-shell photon, for B — > K£ + £~ 
and B — > K*£ + £~ the matrix elements must be known in a wide range of the invariant mass 
squared of the lepton pair: M| + ^_ = [4M|, (Mb — Mk,k*) 2 ]', therefore, the vector meson 
dominance assumption has not negligible consequences on the theoretical outcome. 

An approach based on general features of QCD that allows us to compute the hadronic 
matrix elements in a range of Mf +t _ is provided by three-point function QCD sum rules |TU| 



This method, first employed to compute the pion form factor fill, has been widely applied 



to heavy meson semileptonic decays: for example, in the case of B — > D, D* semileptonic 
transitions, it has been used to compute the Isgur-Wise universal function £(y) and the 



heavy quark mass corrections [12|. Moreover, the decays B — > D**£v, where D** are positive 
parity (cq) meson states, have been analyzed both for finite heavy quark masses [[13] and 
in the limit itlq — > 00, with the calculation of the universal functions Ti(y) and rz{y) 
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analogous to the Isgur-Wise function 

as D(B) — > Tr(p)£v, the various matrix elements have also been computed 
case of B 



For the heavy-to-light meson transitions, such 

in the 



15,16 



iC*7, this approach, employed in [0-1^], has provided us with the prediction 



R = B(B -> K*>y)/B(b -> 57) = 0.17±0.05 0, that agrees with the central value obtained 
from the experimental data in eqs.( |l.l[ )-( |1.2| ). 

In this paper we want to apply the three-point function QCD sum rule method to compute 
the hadronic quantities appearing in the calculation of B — > K^*H + £~ . We shall observe that 
the various form factors parametrizing the relevant matrix elements have common features 
with other heavy-to-light meson transitions, a behaviour whose origin is worth investigating 
in detail pD| . We shall also compare the computed hadronic quantities to the findings of 
lattice QCD, even though these last results are obtained after extrapolations in the heavy 
quark mass and in the momentum transfer. Finally, we shall apply our results to predict 
the invariant mass distribution of the lepton pair in the decays B -> K { *H + l~ and the 
forward-backward asymmetry for B — > K*£ + £~ in the Standard Model. 

The work is organized as follows: in Sec. II we write down the (SM) effective Hamil- 
tonian for the transition b — > s£ + £~, and resume the available information on the Wilson 
coefficients. In Sec. Ill we compute by three-point function QCD sum rules the relevant 
hadronic quantities for B — > K£ + £~; the same calculation is carried out for B — > K*£ + £~ in 



Sec. IV. In Sect. V we study the relations derived by Isgur and Wise [21] and Burdman and 
Donoghue [^2| between rare and semileptonic form factors. Such relations can be worked 
out in the infinite heavy quark mass limit my, — > 00, in the region of maximum momentum 
transfer t; a relevant problem is whether they are satisfied also in the low t region, as it has 
been argued by several authors. We investigate this hypothesis and comment on the role 
of the heavy mass corrections. In Sec. VI and VII we study the transitions B — > K£ + £~ 
and B — > K*£ + £~, respectively. Finally, in Sec. VIII we draw our conclusions. Details 
concerning the calculations are reported in the Appendix. 



II. EFFECTIVE HAMILTON! AN 

The effective AB = — 1, AS = 1 Hamiltonian governing in the Standard Model the rare 



transition b — > s£ + £ can be written in terms of a set of local operators |23] : 

ri 10 

H w = 4 -^V tb V* J2 CiWOiQi) (2.1) 
V2 i=1 

where Gp is the Fermi constant and are elements of the Cabibbo-Kobayashi-Maskawa 



VubV* 



V tb V t * s 



is of the 



mixing matrix; we neglect terms proportional to V u bV* s since the ratio 
order 10~ 2 . The operators Oi, written in terms of quark and gluon fields, read as follows: 

0\ = (s La -f >l b La )(c L i3'y fl c L i3) 

2 = {sLal^bL^iCLplnCLa) 
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A = (sLaYbL^KuLpJuULa) + ... + {b L p1pb La )\ 



5 = {s La ^b La )[{u R ^^u R p) + ... + {b R ^^b R p)} 

O e = (s LQ 7 M 6 L/3 )[(M i?/3 7 M M J? , a ) + ... + (bjip'j^b^)} 
g 

°7 = T7^~o m b{sL a <7^ V b Ra ) Fp V 



16tt 2 



0s= * mb {s L ^£) J 



167T 2 

9 = -^(s La Yb La )h^ 

167T 



C 2 



Ova = -r^(s La Yb La ) h^ld (2.2) 

(a, (3 are colour indices, b RL = — 7T~b, and = -[t' 1 , 7"]); e and g s are the electro- 

magnetic and the strong coupling constant, respectively, and G* in O-j and 0% denote 
the electromagnetic and the gluonic field strength tensor. 0\ and O2 are current-current 
operators, 3 ,...,0 6 are usually named QCD penguin operators, 7 (inducing the radia- 
tive b — ► 57 decay) and 0$ are magnetic penguin operators, O9 and 0\q are semileptonic 
electroweak penguin operators. 

The Wilson coefficients have been partially computed at the next-to- leading order 

in QCD by several groups [E3|-|3E|]. As discussed in ref. in the analysis of B — > X s £ + £~ 
at the next-to-leading logarithmic corrections must be consistently included only in the 
coefficient Cg, since at the leading approximation Og is the only operator responsible of the 
transition b — > s £ + £~ . The contribution of the other operators (excluding 0§ that, however, 
is not involved in the processes we are studying) appears only at the next-to-leading order, 
and therefore their Wilson coefficients must be evaluated at the leading approximation. 
Following |25|] we use in our phenomenological analysis of the decays B — > K^£ + t~ (within 

the Standard Model) the numerical values of the Wilson coefficients collected in Table I. We 

(5) 

choose the scale fi = 5 GeV ~ m&, = 225 MeV and the top quark mass m t = 174 GeV 
from the CDF measurement ||271| . The coefficient Cg, which is evaluated at the next-to- 



leading order approximation, displays a dependence on the regularization scheme, as it can 
be observed in Table I comparing the result obtained using the 't Hooft-Veltman (HV) and 
the Naive Dimensional Regularization (NDR) scheme. Such dependence must disappear 
in the decay amplitude if all corrections are taken into account. We shall include in our 
analysis the uncertainty on Cg as a part of the theoretical error. In Table I it can also be 
observed that the coefficients of O3 — 0$ are small (0(1O~ 2 )); therefore, the contribution of 
such operators can be neglected, and the analysis can be carried out considering only the 
operators Oi, 2 and to 7 , Og and O w . 

The various extensions of the Standard Model, such as models involving supersymmetry, 
multiHiggs and left-right models, induce two kind of changes in the low energy Hamiltonian 
(|2.1|) : first, the values of the coefficients Cj are modified as an effect of additional virtual 
particles in the loop diagrams describing the b —>■ s transition, and, second, new operators 
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can appear in the operator basis, such as operators with different chirality of the quark 

e e 2 - 

current with respect to 7 - O 10 , e.g., 0' 7 = -^^ , m b (s R a flu b L )F flI/ , 0' 9 = j^(s R 'y' J, b R )£'y li £ 

e 2 

and O[ = —^(s R Yb R )h^ b L 

This rich structure justifies the interest for B — > K^*H + £~ , where operators of different 
origin act coherently in determining rates, spectra and asymmetries. For example, it could 
be interesting to search for the effects of possible interactions that produce a coefficient 
Cj with opposite sign |5|,[7j. In this work we shall not analyze such new effects, limiting 
ourselves to studying the above processes within the theoretical framework provided us by 
the Standard Model. However, it is worth stressing that our results for the hadronic matrix 
elements of the operators appearing in ( [2.1] ) represent a complete set of quantities also for 
the analysis of the decays B — ► K^£ + £~ in a context different from the Standard Model. 



III. FORM FACTORS OF THE DECAY B -»■ Ki + £~ 

The matrix elements of the operators Oi, O2 and O7, O9 and O\o in eq.( [2.2| ) between 
the external states B and K can be parametrized in terms of form factors as follows: 

< K(p')\s^b\B(p) >= ip + v)^ 2 ) + B 2 K Q» (Fo(q 2 ) ~ F^q 2 )) (3.1) 



q 2 



q = p-p', Fi(0) = F (0)) and 



<K(p')\sia, u q»b\B(p) >= [(p + p\q 2 - (M| - M|)g M ] • ( 3 " 2 ) 

The heavy-to-light meson form factors i*\ and Fo appear in the calculation of two-body 
nonleptonic F — > FX decays, if the factorization approximation is adopted; neglecting 
5F(3)^ breaking effects, they govern the semileptonic decay F — > 7r£z/. 

Fi and F have already been studied by three-point QCD sum rules [|To],EHJ. In the 



following we describe in detail the calculation of Ft; for the sake of completeness, we also 
report the results for Fi(q 2 ) and F (g 2 ) using a unique set of parameters and adopting a 
coherent numerical procedure, in order to have at our disposal a consistent set of form 
factors. 

To compute Ft within the QCD sum rule approach we consider the three-point correlator 





U a ^(p,p',q) = i 2 j dxdy e^'-v-^) < |T[J*(y)J^(0)jf (z)]|0 > (3.3) 



of the flavour changing quark current = sia^b and of two currents J^iv) an d Ji{%) 
with the K and B quantum numbers, respectively: J^iv) — Q.{y)lal5 s {y) an d Ji{%) = 
6(x)i7 5 g(x). The correlator H afll/ can be expanded in a set of independent Lorentz structures: 

= ip'MpI - iw<jn + i £ <l n(n) (3.4) 

n 
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where II and IT") are functions of p 2 , p' 2 and q 2 , and a^ v are other tensors set up using the 
vectors p and p' and the metric tensor g^ v . 

Let us consider II. To incorporate the quark-hadron duality, on which the QCD sum 
rule approach is based, we write down for U(p 2 , p' 2 , q 2 ) a dispersive representation: 

]_ r+oo r+oo p(s,s',q 2 ) 

Il(p 2 ,p' 2 ,q 2 ) — — / ds ds' - ^-r 2 - — + subtractions (3.5) 

7T 2 Jml J ml (s - p 2 )(s' - p' 2 ) 

in the variables p 2 and p' 2 corresponding to the B and K channel, respectively. In the region 
of low values of s, s' the physical spectral density p(s, s', q 2 ) contains a double 5-function 
term corresponding to the transition B — > K, and therefore the function II can be written 



as 



{M 2 -p 2 ){M 2 K -p> 2 ) + <k 2 Jd {s-p 2 ){s' -p> 2 ^ > " >l 



where the residue R is given in terms of the form factor i^(g 2 ) and of the leptonic 
constants fx and defined by the matrix elements < 0|g7 M 7 5 s|i^(p') >= ificPu 
and < 0\qi>y 5 b\B(p) >= f B M 2 B /m b (we put m q = 0): R — HF T (q 2 ) with H = 
—2fKfBM^/mb(MB + Mk)- The integration domain D in (|3.6|) , where higher resonances 
with the same B and K quantum numbers contribute to the spectral density p, starts from 
two effective thresholds s and s' Q . 

Also the perturbative contribution to EE, computed for p 2 — > — oo and p' 2 — ► — oo, can 
be written as eq. ( |3.5| ). Moreover, considering the first power corrections of the Operator 
Product Expansion of the correlator (|3.3|) we get the following representation: 



(3.7) 



(s, s', q 2 ) is the perturbative spectral function; the two other terms in Q3.7|) , expressed as 
a combination of vacuum expectation values of quark and gluon gauge-invariant operators of 

dimension 3 and 5, respectively: < qq > and < qoGq >=< gsqa^G^—q >, parametrize 

the lowest order power corrections. The expressions for pQ° D and d 5 can be found in 
Appendix A, eqs. (|A2|) - (|A4|) ; in this particular case d^ vanishes. 

We now invoke the quark-hadron duality, i.e. we assume that the physical and the 
perturbative spectral densities are dual to each other, giving the same result when integrated 
over an appropriate interval. Assuming duality in the region D of the hadronic continuum 

/ dsds'{p had (s, s', q 2 ) - p QCD (s, s', q 2 )} = (3.8) 

JD 

we derive the sum rule for Ft' 

HF T (q 2 ) 1 f , , , p^ CD (s,s',q 2 ) 

{Ml-p 2 ){M 2 K -p< 2 ) =- 2 L dsdS (s-p 2 )(s>-p> 2 ) +d * < ^ > +d * < qaGq > + - ' 

(3.9) 
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where D' is the region corresponding to the low-lying B and K states: m 2 < s < sq, s'_(s) < 

- mi 



s' < s',(s) with s'i(s) = m 2 + 



and 



?±yo) - ,,o s -r ^ 2 {m\ + m\ - q 2 ) ± yj (m 2 b + m 2 s - q 2 ) 2 - Am 2 b m 2 s 

s+ < s o- The effective thresholds So and s' can be fixed from the QCD sum rule analysis of 
two-point functions in the b and s channels. We get s from the calculation of and s' 
from the expected mass of the first radial excitation of the kaon. 

An improvement of the expression in fl3.9| ) can be obtained by applying to the left and 
right hand sides the SVZ-Borel transform, defined by 

I l e -m 2 /M 2 

EM \m 2 - V 2 Y = (n- 1)! (M 2 )™ ' 

both in the variables — p 2 and — p' 2 ; M 2 is a new (Borel) parameter. This operation has 
the advantage that the convergence of the power series is improved by factorials; more- 
over, for low values of M 2 and M' 2 the possible contribution of higher states in eq.( ft.9|) is 
exponentially suppressed. The resulting Borel transformed sum rule for Ft reads 

HF T (q 2 )e- M B/ M2 - M y M ' 2 = 1 / dsds'p QCD (s, s', tfy-'/M*-'/**" 

7T 2 JD' 

+ \d 3 <qq> +d 5 < qoGq > } e -</M 2 -m*/M>* _ ^ n ^ 



From eq. (|3.11|) the form factor ^(g 2 ) can be derived, once the value of the Borel pa- 



rameters M 2 and M' 2 is fixed. This can be done observing that, since M 2 and M' 2 are 
unphysical quantities, Ft must be independent on them (stability region of the sum rule); 
moreover, the values of M 2 and M' 2 should allow a hierarchical structure in the series of the 
power correction, and a suppression of the contribution of the continuum in the hadronic 
side of the sum rule. 

In our numerical analysis we use the values for the quark condensates (at a renormaliza- 
tion scale /i ~ 1 GeV) |Tl| : 

< qq > = (-230 MeVf 

< q9s^ u G%^-q > =m\ <qq> (3.12) 

with m 2 , = 0.8 GeV 2 . Notice that the numerical results do not change sensitively if the 
condensates are evaluated at higher scales using the leading-log approximation for their 
anomalous dimension. 

As for the quark masses and leptonic constants, we use: m s = 0.175 GeV, m^ = 4.6 GeV, 
f K = 0.16 GeV and f B = 0.18 GeV. The thresholds s and s' are chosen in the range: 
so = (33 — 36) GeV 2 and s' = (1.4 — 1.6) GeV 2 , with the Borel parameters kept fixed to 
the values M 2 = 8 GeV 2 and M' 2 = 2 GeV 2 . 

Putting these parameters in eq. ( |3.11| ) we obtain the form factor Ft depicted in fig.l, 
where the different curves correspond to different choices of the thresholds so and s' Q . In the 
sum rule, the perturbative term is a factor of 4—5 times larger than the D = 5 contribution, 
and the integral of the spectral function over the region D' gives more than 60% of the 
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result of the integration over the whole region of the dispersion relation Q3.5|) . The duality 
window, where the results become independent of the Borel parameters M 2 and M' 2 , starts 
at M 2 ~ 7 GeV 2 and M' 2 ~ 1.7 GW 2 ; varying M 2 in the range 7-9 GeV 2 and M' 2 in the 
range 1.7 — 2.5 GeV 2 the results change within the bounds provided by the different curves 
depicted in fig.l. 

The same analysis can be applied to the form factors F\ and Fq using the flavour- 
changing vector current J M = s^^b in the correlator (|3.3j) and studying the projection gTI a(U 



to derive Fq. We report in Appendix A the relevant quantities appearing in the sum rules 
for F 1 and F ; the difference with respect to [T5I] , as far as F 1 is concerned, is that we 
keep all terms proportional to powers of the strange quark mass m s . In the calculation of 
both the form factors, the contribution of the perturbative term and of the D = 3 term have 
comparable size, whereas the D = 5 term is one order of magnitude smaller; the contribution 
of the resonance in the hadronic side of the rule is nearly equal to the contribution of the 
continuum. We obtain the form factors Fi(q 2 ) and F Q (q 2 ) depicted in fig.l. Also in this case 
the Borel parameters can be varied in the range M 2 = 7—9 GeV 2 and M' 2 = 1.7 — 2.5 GeV 2 ; 
the results change within the region corresponding to the different curves depicted in fig.l 
for each form factor. 

We observe a different q 2 dependence for the various form factors. In the range of q 2 we 
are considering (0 < q 2 < 13 — 15 GeV 2 ) Fi follows a simple pole formula: 

F(q 2 ) = (3.13) 
1_ M 2 



with Fi(0) = 0.25 ± 0.03 and M Pl ~ 5 GeV. A fit to the formula (|3TT3| ) for F gives 



the result Mp ~ 7 GeV. The same formula, applied to Ft would give Ft — —0.14 and 
Mp ~ 4.5 GeV. Therefore, only the dependence of the form factor Fi(q 2 ) does not contradict 
the polar behaviour dominated by B*, which is the nearest singularity in the t— channel, as 
we would expect by invoking the vector meson dominance (VMD) ansatz. The form factor 
F increases softly with q 2 and, as already observed in ||28|| , the fitted mass of the pole is 
larger than the expected mass of the physical singularity, in this case the J p = + bs state. 
As for Ft, the VMD ansatz would predict a polar dependence, with the pole represented 
by B*; on the other hand, we observe that Ft can be related to F\ and F by an identity 
obtained by the equation of motion: 

F T (q 2 ) = (M B + M K ){m b + m / ^"/ 1 ^ 1 ; (3.14) 

q z 



eq. (|3.14|) is in agreement with the computed form factor Ft displayed in fig. 1, and therefore 



we can use the double pole model: 



Ml') = ? x (3-15) 



1 




with Ft(0) = —0.14 ± 0.03 and Mp 1 and Mp given by the fitted values of the mass of the 
poles of Fx and F , respectively. 



S 



It is interesting to observe that information on the possible form of the q 2 dependence 
of the form factors can be derived by studying the limit — > oo. In this limit, at the zero 
recoil point where the kaon is at rest in the B meson rest frame, it is straightforward to show 
that the parametric dependence of the form factors on the heavy meson mass Mb is given 
by: Fi(q 2 max ) ~ \JMb and F (q^ nax ) ~ 1 /V Mb [[HJ]. Both these scaling laws are compatible 
with the constraint -Fi(O) = Fq(0) and with a multipolar functional dependence 

F^ 2 ) = ( Fli °Ln, (3.16) 



Ml 



if n\ = no + 1. Thus, in the limit mj — > oo, to a polar form factor Fx(q 2 ) corresponds a 
nearly constant form factor Fo(q 2 ). The outcome of QCD sum rules is in agreement with this 



observation |p9| ; the observed increasing of F would be due to subleading terms contributing 
at finite m^. 

Let us now compare our results with the outcome of different QCD based approaches. 
In the channel B — > 7r the form factor F± has been computed by light-cone sum rules |30| . 



with numerical results in agreement, at finite b-quark mass, with the outcome of three point 
function sum rules. 

As for lattice QCD, both Fi and F have been computed at large q 2 [31], and data show 



that Fq has a flat dependence on the momentum transfer, whereas Fx increases with q 2 . 

The full set of form factors Fx, F Q and Ft by these other methods is still missing; the 
complete comparison of our results with such different approaches could help in understand- 
ing the drawbacks and the advantages of the various methods; this would shed light on the 
issue of decays such as B — > Triu that are of interest as far as the measurement of V u b is 
concerned. 



IV. FORM FACTORS OF B — > K*£+£- 

The form factors parametrizing the hadronic matrix elements of the transition B — > 
K*£ + £~ can also be computed by QCD sum rules by considering a three-point correlator 
with the interpolating current for K* represented by the vector current j£' '(y) = q(y)j a s(y). 
Let us define the B — > K* matrix elements: 

< k*( p ', 6)|^(i - i*)b\B{p) > = e^yy 3 2V{q2) 



— I 



M B + M K , 

el(M B + M K *)Ax{q 2 ) - (e* • q )(p + ^—^S. 



and 



(e* • g)^l(A 3 (g 2 ) - A (q 2 ))q, 



(M B + M K *) 
(4.1) 
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< K*(p\e)\sa^^p^b\B(p) > = ie^ a/3 e*»p a p'P 2 T\(g 2 ) + 



+ 



e ;(M B -M 2 K *)-(e*-q)(p + p%] T 2 (q 2 ) 



+ (e* ■ q) 



g "~ M 2 -M 2 J P + P ^ 



T 3 (q 2 ). (4.2) 



A3 can be written as a linear combination of A-y and A 2 . 

2 M B + M K « 2 ^ Mb-Mk, 2 

with the condition A 3 (0) = A (0). The identity 0^75 = — -e Ati , a/ 3<7 a/3 (e i23 = +1) implies 

that Ti(0) = T 2 (0). 

The form factors T\(q 2 ) and T 2 (g 2 ) can be derived by the correlator 

U afl {p,p',q)=i 2 J dxdye&'-y-^ < 0|T[jf *(y) J M (0) J 5 B (x)]|0 > , (4.4) 

1+75 

with J M = scr^ — - — Expanding Il QAt in Lorentz independent structures 

ria/i = ie Q , Atp/ 3p p p //3 ili + g QAt il 2 + ot/zer structures in p, p' (4.5) 
we get T\ and T 2 from ITi and Il 2 , respectively. The sum rules have the same structure of eqs. 



( p.9[) , (|3.11|) , with the perturbative spectral functions p(s, s', q 2 ) and the power corrections d% 
and c?5 reported in Appendix B. The only difference with respect to the kaon case is the value 
of the K* leptonic constant, defined by the matrix element < 0\qj fl s\K*(p, e) >= fK*M K *e^, 
with f K , = 216 MeV. 

In fig. 2 we depict the form factors Ti(q 2 ) and T 2 (g 2 ) obtained choosing the threshold 
Sq in the range 1.6 — 1.8 GeV 2 and the other parameters as in the previous section. In 
the sum rule for both the form factors the perturbative term does not dominate over the 
non-perturbative ones: at q 2 = it represents 30% of the quark condensate contribution, 
and is nearly equal to the D = 5 term. However, it rapidly increases with the momentum 
transfer, and at q 2 = 15 GeV 2 it is equal to the contribution of the D = 3 term, whereas 
the D = 5 contribution is an order of magnitude smaller. 

Concerning the form factor T3, we observe that it contributes, together with T\ and T 2 , 
to other invariant functions in (|4.5| ) and, in principle, it also could be obtained by a sum 
rule. However, since it can be related to Ax, A 2 and A Q by applying the equation of motion: 

T 3 (q ) = M K *(m b - m, 5 (4.6) 

qz 

we prefer to use this expression to determine it, considering that this procedure is successful 
for F T (q 2 ). 

The form factors V and Ai can be obtained by studying the correlator ( |4.4[ ) with a vector 
J* = an d an axial = 57^756 flavour changing current, considering the projection 



10 



q^J^ to derive A . We collect in Appendix B the complete expressions appearing in the 
relevant sum rules for all the form factors, excluding A , whose expressions can be found in 
1 52| ; also in this case the difference with respect to [[DJ is that we include all powers of the 
strange quark mass. 

Using our set of parameters we get V(q 2 ), Ai(q 2 ), A 2 (q 2 ) and A (q 2 ) depicted in fig. 3, 
and, using (|4.6|) , the form factor T3 in fig. 2. 

As it happens for T\ and T 2 , also in the sum rules for V , A\ and A2 the perturbative 
term, at q 2 = 0, is smaller than the D = 3 contribution; the relative weights of the various 
contributions change with the momentum transfer, and at q 2 = 15 GeV 2 the D = and 
D = 3 terms have comparable size. As it happens for the B —>■ K form factors, the chosen 
values of M 2 and M' 2 , M 2 = 8 GeV 2 and M' 2 = 2 GeV 2 , are within the duality window 
where the results are independent of the Borel parameters. Also in this case, varying M 2 
and M' 2 in the ranges M 2 = 7 — 9 GeV 2 and M' 2 = 1.7-2.5 GeV 2 , the final results change 
within the same uncertainty coming from the variation of the continuum threshold. 

Considering the results displayed in figs. 2 and 3, we collect the form factors T iy V and Ai 
in three sets, according to their functional dependence on the momentum transfer. In the first 
set we include 7\, V and Aq, that display a sharp increasing with q 2 . It is possible to fit them 
with a polar q 2 dependence eq.( |3.13| ) (as observed also in (|Tj||3^]) with: T^O) = 0.19 ± 0.03 
and M P ~ 5.3 GeV, V(0) = 0.47 ± 0.03 and M P ~ 5 GeV, 4,(0) = 0.30 ± 0.03 and 
M P ~ 4.8 GeV (the difference with respect to the value 7\(0) = 0.17 ± 0.03 in ref. is 
due to the effect of the strange quark mass, that here has been included). 

The error on the mass of the pole is correlated to the error on the form factor at q 2 = 0, 
and it can be estimated of the order of 200 — 300 MeV. The relevant result is that the masses 
of the poles are not far from the values expected by the dominance of the nearest singularity 
in the t— channel: Mp = Mp* for T\ and V , Mp = Mp g for A . We stress that the fit is 
performed in a range of values of q 2 where the QCD calculation can be meaningfully carried 
out, therefore large momentum transferred [q 2 > 15 GeV 2 ] are not taken into account. 

In the second set of form factors we include T 2 , T3 and A\. They softly decrease with 
q 2 : Fi(q 2 ) = ^(0)(1 + f3q 2 ), with T 2 (0) = T\(0) and (3 = -0.02 GeV' 2 , T 3 (0) = -0.7 and 
P = 0.005 GeV~ 2 , A x (0) = 0.37 ± 0.03 and (3 = -0.023 GeV~ 2 with the error on (3 at the 
level of 10%. The dependence of T 3 is related to A\, A 2 and A . 

The last form factor, A2, linearly increases with q 2 : ^(0) = 0.40 ± 0.03 and (3 = 
0.034 GeV~ 2 . A fit to a polar dependence for this form factor would give Mp > 7 GeV for 
the mass of the pole. 

The parameters of all the form factors are collected in Table II. Albeit the form factors 
have been computed in a well defined range of momentum transfer, once their functional 
q 2 dependence has been fitted and the parameters determined, we extrapolate them up to 
q^ax- This procedure cannot be avoided within the method of QCD sum rules, where large 
positive values of q 2 are not accessible since there is a region where the distance between 
the points x, y and in the correlator, which is the initial ingredient of this approach, is 
large, and therefore the standard OPE cannot be used; this is shown by the occurrence of 
singularities in the correlator when q 2 is close to q^^. 

As for the computed dependence on the momentum transfer, is worth reminding that 
deviations from the VMD expectations for the form factors A\ and A2 have been already 
observed in the literature, first in the D — > K*£u JT5J channel and then for B — > plv 
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|16|| . Here we find a kind of common feature, i.e. all form factors deviating from the polar 
dependence (excluding Ft) seem to depend linearly on the momentum transfer, with small 
(positive or negative) slopes. 

It is interesting that also for Ti(q 2 ) and T 2 (q 2 ) we can use the argument developed in 
the previous section concerning the limit vrt}, — > oo: since Tx(q^ nax ) ~ \JMb and T 2 (q 2 nax ) ~ 

\/\J Mb, the constraint Ti(0) = T 2 (0) can be fulfilled by a multipolar q 2 dependence if 
rii = n,2 + 1 in eq. ( |3.16| ). 

At zero momentum transfer our results numerically agree with those obtained by the 
method of light-cone sum rules |£J , within the errors and taking into account the different 
choices of the input parameters. In [33] it has also been observed that T 1; V and A\ have 
different functional dependencies on q 2 ; the difference with respect to our case is that the 
slopes are larger than those obtained from three-point sum rules; in particular, the form 
factor A\ increases with q 2 . The origin of this discrepancy should be investigated. 



The form factors T\ and T 2 have been computed by lattice QCD [|34| , |3~5"f near the point 
of zero recoil and for the mass of the heavy quark smaller than rrib, due to the finite size 
of the available lattices; therefore, the results at q 2 = and for a realistic value of m;, are 
obtained after an extrapolation in the momentum transfer and in the heavy quark mass. 
Also in this case, in the region of large values of q 2 , the form factor 7\ increases rapidly 
with the momentum transfer, whereas T 2 is quite flat. As for the analytic q 2 behaviour 
obtained from lattice calculations, it seems to us that larger lattices are needed to enlarge 
the range of momentum transfer where the measurements can be performed, in order to 
clearly disentangle different possible dependencies of T\ and T 2 (e.g., dipole versus pole or 
pole versus constant). 



V. RELATIONS BETWEEN RARE AND SEMILEPTONIC B DECAY FORM 

FACTORS 



In the limit — > 00 Isgur and Wise |21|] and Burdman and Donoghue |22] have derived 
exact relations between the form factors Ft, Ti in eqs.(EO), ( ^.2| ) and the form factors i^, 
V, Ai in eqs. (|3.1| ), (f4.1|) . These relations can be easily worked out observing that, in the 
effective theory where the 6-quark mass is taken to the infinity, the equation 7°6 = b is 
fulfilled in the rest frame of the B meson. 

In our parametrization such relations can be written as follows, near the point of zero 
recoil (q 2 ~ q 2 max = (M B - M u 



F T (q 2 



M B + M K 



2M 



B 



^i(9 2 )-(Ml-Mi; 



F (q 2 )-F 1 (q 2 



(5.1) 



Uq 2 



M B + M K . 



AM, 



B 



+ 



Ml - Mg + q 2 
4M B (M B + M K * 



-V{q 2 



(5.2) 



T 2 (g 2 



(Ml - M 2 K , + q 2 
4Mb (Mb ~ M K , 



-Mq 2 ) + 



X(M 2 ,Mh,q 2 



AM B (M B - M K *)(M B + M K * 



-,V{q 2 



(5.3) 
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T 3 tf) = Ml + 3Mf. - g 2 + M K * Mq2) + M K *(M 2 - M%.) A 3 (q 2 ) - A (q 



4M B (M B + M K *) ^ J 2M B 2M b q 2 

(5.4) 

where A is the triangular function. 

It has been argued by several authors that the relations fl5.1| )-(|5Tl) could also be valid 



at low values of q 2 |p2]| , although a general proof has not been found in support of this 
hypothesis. 

Using the form factors computed by QCD sum rules in the previous Sections, it is possible 
to check eqs.( [5~Tl) -( |5.4|) . In fig.4 we plot the ratio 1Z = F/F IW in the case of Ft, Ti, and 
T2, as a function of q 2 , in the range of momentum transfer where the calculation has been 
carried out. We observe that the relations between the various form factors are verified at 
different level of accuracy. 

In the case of Ft the ratio TZ differs from unity at the level of 25 — 30%, including the 
uncertainty coming from the errors of the various parameters. In particular, at q 2 = we 
have Ft/F t w = 0.7 ± 0.1. The situation is different for the ratios concerning T\ and T2, 
that differ from unity at the level of 10 - 20%: at q 2 = we have Ti/T/ W = 0.94 ± 0.05 and 
T 2 /T 2 /w = 1.12 ±0.05. 

These results support the argument put forward in [17| on the validity of the Isgur-Wise 



relations, in the limit m& — > 00 also at small values of q 2 ; they also can be well compared 



to the outcome of light-cone sum rules, obtained for Ti at a finite m& ||33|| . The conclusion 
is that the b quark is near to the mass shell also when the recoil of the light hadron is large 
with respect to mj,, with I/777.& corrections that do not appear to overwhelm the effect. 

The relations (|5.1|)-(]5T4|) could be used to perform a model independent analysis of the 
decays B — > K^*'£ + £~ employing experimental information (when available) on the form 
factors of the semileptonic transition B — > ptv ||36|| . In particular, since (|5.1| )- (|5.4| ) are valid 
on general grounds in the large q 2 region, it has been proposed to perform the analysis in 
the range of large invariant mass of the lepton pair, e.g. Mp+p- > 4 GeV. 

Albeit in principle correct, we feel that, from the experimental point of view, the proce- 
dure of extracting the semileptonic B —>■ p form factors near zero recoil will be rather difficult, 
with large uncertainties in the final result. The problem is not avoided by the possible choice 
of using the form factors of the semileptonic transition D — > K*£ + u, and then rescaling them 

V Bp (q 2 ) /Mb 

according to the their leading dependence on the heavy mass, i.e. , rr .^, m ^ x x = \ 

V DK WmaJ V M D 

etc. (neglecting SU(3)p and a s corrections). As a matter of fact, in such procedure the 
next-to-leading mass corrections could be large and not under control. Finally, as we shall 
see in the next section, the differential branching ratios of B — > K^*'£ + £~ at large q 2 are 
small, and therefore the experimental errors are expected to be sizeable. For this reason we 
prefer to propose an analysis of the decay extended to the full range of q 2 , using hadronic 
quantities determined in a well defined theoretical framework. The dependence on the com- 
putational scheme will be reduced once the different form factors have been computed by 
different QCD calculations, and the whole information collected in a unique set of form 
factors. 
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VI. DECAY B -> K£+£~ 



We can now compute the invariant mass squared distribution of the lepton pair in the 

decay B -»■ KPtr: 



1536tt 5 
F T {q 2 ) x 



M B + M K 



+ 



CioF^g 2 ) 



Ml, 



+ 



MS, 



Ml, 



Ml, 



3/2 



(6.1) 



(g = M £+£ _). The contribution of the operators O7, Og and O10 is taken into ac- 
count in the terms proportional to C7, Cg and C10. The operators 0\ and O2 provide a 
short distance contribution, with a loop of charm quarks described by the function h(x, s) 
(x = m c /m b , s = q 2 /mf*j 



h(x, s) = — 



4 2 8 16 x 2 4 /4x 2 
- nx - 27 - y— + gV — 



/ 4x 2 \ /4x 2 
1 ( 2 + — 1 arctg I — 



-1/2 



(6.2) 



if s < 4a; 2 and 



/i(x, s) 



4, 9 8 16x 2 2 / 4x 2 / 4x 2 

-Inx 2 1 — a/1 2 H 

9 27 9s9V s s 



in 



1 + Jl -4x 2 /s 



1 - Jl -4x 2 /s 



(6.3) 



if s > 4x 2 ; the imaginary part in ( |6.3D comes from on-shell charm quarks. 0\ and O2 also 
provide a long distance contribution, related to cc bound states (J/ip^ip 1 ) converting into 
the lepton pair £ + £~~ [^,^]. This contribution can be described in terms of the J/i/j and 



ip' leptonic decay constants < 0|c7 M c|^i(e, q) >= e^f^M^ and of the full J/ip and ip' decay 
widths T^. We derive fy. from the experimental branching ratio ipi — > £ + £~; in this way 
the whole contribution of 0\ and O2 can be taken into account by modifying the coefficient 
CqintoC Q e// : 



C 9 + (3d + C 2 



h{x lS ) + kj: 7:T ^ 



l+l-)M^ 



(6.4) 



Iqi-Ml + iM^ 
Kipi transition is computed by factorization, the parameter k is 



If the nonleptonic B — 

given by k = — — cs ; the sign between the short distance and the long distance term 



a 2 \V t * s V th \ 

in ( |6.4|) can be fixed according to the analyses in ref. 
appropriately chosen in order to reproduce the quantity: 



In ref. [f| the value of k is 
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B(B -> K£ + £~ ) = V -> ^iK)B(^i -> £+r) ~ 7 10~ 5 [39]. (6.5) 

res — ' 

i=l 

3 

This can be done by choosing k ~ (1.5 4- 2) x — . Notice that, since the J lip and ip' 

cr 

resonances are narrow, their contribution modifies the dilepton spectrum only in the region 
close to M 2 +l . = Mj /f , M|,. 

As input parameters we choose the ratio m c /mi, = 0.27 — 0.29 and the value of the CKM 
matrix element |Vt s | ~ 0.04; a different value for \Vt s \ only modifies the prediction of the 



branching ratio, leaving unchanged the shape of the spectrum |40 



We depict in fig. 5 the obtained invariant mass squared distribution of the lepton pair 
in B — > K£ + £~~ . In the same figure we also plot the spectrum obtained considering only 
the short distance contribution, that gives the branching ratio (using tb = 1.5 10~ 12 sec 
for the B- meson lifetime) B(B -> K£ + £~)\ sd ~3x 10~ 7 |V ts /0.04| 2 , to be compared to 
the experimental upper limit (obtained excluding the region near Jjip and ip r ) B(B~ — > 
K~fi + [i~) < 0.9 x 10~ 5 (at 90% CL) PIj^2| . The uncertainty coming from the two possible 
values of Cg in Table I is less then 1% and does not have relevant consequences on the 
predicted branching ratio and on the invariant mass distribution. 

From the experimental point of view, the measurement of the spectrum in fig. 5 is a non 
trivial task; hopefully, it will be possible to obtain experimental results from the future 
dedicated e + e~ colliders. The important point to be stressed is that, in the distribution 
depicted in fig. 5 the theoretical uncertainty connected to the hadronic matrix element is 
reduced to a well defined QCD computational scheme (QCD sum rules), so that in the 
studies of the effects of interactions beyond the Standard Model the hadronic uncertainty 
plays no more a major role. 



VII. DECAY B -> K*£+£- 

A great deal of information can be obtained from the channel B — > K*£ + £~ investigat- 
ing, together with the lepton invariant mass distribution, also the forward-backward (FB) 
asymmetry in the dilepton angular distribution; this may reveal effects beyond the Standard 
Model that could not be observed in the analysis of the decay rate. 

A FB asymmetry in the dilepton angular distribution is a hint on parity violation. Since 
the decay B — > K*£ + £~ proceeds through 7, Z and W intermediate bosons, we expect a 
different behaviour in the various q 2 kinematical regions. In the region of low q 2 , the photon 
exchange dominates, leading to a substantially vector-like parity-conserving interaction; as 
a consequence, we expect a small asymmetry. On the other hand, when q 2 is large, the 
contribution of Z and W exchange diagrams becomes important, and the interaction acquires 
the V-A parity violating structure, leading to a large asymmetry. As already observed in 
ref. [^] this pattern strongly depends on the value of the top quark mass, and the penguin 



diagrams with Z exchange and the W box diagram are expected to overwhelm the photon 
penguin diagram in correspondence to the measured mj. Moreover, since the FB asymmetry 
is sensitive not only to the magnitude of the Wilson coefficients, but also to their sign |J, it 
can be used to probe the values predicted by the Standard Model. 

Let us define 8g as the angle between the £ + direction and the B direction in the rest 
frame of the lepton pair. Since, in the case of massless leptons, as we assume, the amplitude 
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can be written as sum of non interfering helicity amplitudes, the double differential decay 
rate reads as follows: 

d 2 T Gl\V tb V t :\W\^(MlMh, q 2 ) r 2 

= 913^6 m| x t«» Ml+ 



dq 2 dcos0£ 



+ q 2 [(1 + cos^) 2 (^ + A K ) + (1 - cosO t )\A L _ + A*)] } (7.1) 

where Al corresponds to a longitudinally polarized K*, while A + \{ represent the con- 
tribution from left (right) leptons and from K* with transverse polarization: e± = 



We obtain 
1 



A L = 



M 2 K , 



Bi(Mg — M\, — q 2 ) + B 2 \ + D 1 (M 2 B - M 2 K , - q 2 ) + D 2 \ (7.2) 



and 



A L ± = \\ 1 ' 2 {A-C)T{B l -D l ) 



(7.3) 



^ = |A 1 / 2 (A + C7)=f(S 1 + £> 1 )| 2 , 



(7.4) 



where A = A(M|, M|-*, g 2 ). The terms A,C,Bi,Di contain the short distance coefficients, as 
well as the form factors: 



A = ^Am b T 1 (q 2 )+C 9 



V{q 2 



M B + M K * 



(7.5) 



C = C 
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V(q 2 



M B + M K * 



(7.6) 



B 1 = ^4m b T 2 (q 2 )(M 2 - M%.) + C 9 A l {q 2 ){M B + M K *) 



(7.7) 



5, = - 



C< 7 

q 



2 4m b {Uq 2 )+q 2 



+ C 9 



{M% -M 2 K *)) M B + Mx 



D 1 = C 10 A 1 (q 2 ) (M B + M K * 



Do = -C 
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M B + M K * 



(7.8) 
(7.9) 
(7.10) 



The FB asymmetry is defined as 
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Jo 



d 2 T f° d 2 T 

-dcosQo — \ - n , —dcosO/ 



aFB ( 2\ _ Jo dq 2 dcos6e J-i dq 2 dcosQi , 711 , 



d 2 r J n f o d 2 r 1 r 

-dcos&e + / - - —dcostli 



o dq 2 dcos0£ J-i dq 2 dcos6. 



thus we have 



{q ' 4A L + 2q 2 (AL + A* + AL + A R ) ' K ' 

A FB (q 2 ) is depicted in fig. 6; it is consistent with the prediction of low asymmetry in the 
small q 2 region and high asymmetry for large q 2 . The analysis of the individual shapes of 
the helicity amplitudes (neglecting the long distance contribution) shows that A+ and 
have comparable size, and therefore there is a cancellation of their contribution in eq. ( [iM2j ); 
moreover, they are small with respect to A_' R . In the region of large M 2 +e -, A 1 ^ dominates 
over A R , whereas the situation is reversed for low dilepton invariant mass squared, and this 
is the reason of the small positive asymmetry appearing in fig. 6 for Mf +E - < 3 GeV 2 . It is 
interesting to observe that such positive asymmetry depends on C 7 , and that it disappears 
if C 7 has a reversed sign. 

The invariant mass squared distribution of the lepton pair is depicted in fig. 7, where 
the short distance contribution is separately displayed. The predicted branching ratio is 
B(B — ► K*£ + £~)\ sc i = 1 x 10 -6 | V^/0.04| 2 , to be compared to the experimental upper limit: 
B(B° -»■ K*°fi + {i-) < 3.1 x 10- 5 (CLEO II) and B(B° -> i^V+zO < 2.3 x lO" 5 (UA1) 
(at 90% CL) obtained excluding the region of the resonances J /if) and if)' ||4T| , f44|| , f39fl . Also 
in this case the uncertainty on Cg does not have relevant consequences. 

The interesting observation is that, for low values of the invariant mass squared, the 
distribution is still sizeable, an effect that could be revealed at future S-factories such as 
the Pep-II asymmetric e + e~ collider at SLAC. 



VIII. CONCLUSIONS 

In this paper we have analyzed some features of the rare decays B — > K£ + £~ and 
B — ► K*£ + £~ within the theoretical framework provided by the Standard Model, using an 
approach based on three point function QCD sum rules to compute the relevant hadronic 
matrix elements. 

Albeit QCD sum rules have their own limitations (finite number of terms in the Operator 
Product Expansion of the correlators, values of the condensates, validity of the local duality 
assumption), we believe that the obtained results are meaningful from the quantitative point 
of view. 

There is a quite good agreement with independent QCD methods (lattice QCD, light- 
cone sum rules) for few quantities computed by the various approaches. The calculations 
of the remaining quantities (F , T i} A ) by the other two methods is required in order to 
complete the overview on the various results. 

We have used our results to test some relations among the computed form factors which 
hold in the infinite heavy quark limit, but that are expected to hold also for low values of 



17 



q 2 and for finite b mass. We have found that the different form factors satisfy with different 
accuracies these relations, which can be explained by a different role of the 1 / corrections. 

As for the decays we have analyzed in the present paper, within the Standard Model they 
are expected with branching ratios of the order 10~ 7 (B — > K£ + £~) and 10~ 6 (B — > K*£ + £~), 
with peculiar shapes of the invariant mass of the lepton pair and of the FB asymmetry. Any 
deviation from the above expectations would be interpreted as a signal of deviation from 
the Standard Model. Interesting experimental data are therefore expected from current and 
future e + e~ colliders in this exciting sector of the heavy flavour physics. 
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APPENDIX A: B i^+r 



The three-point function QCD sum rule for the form factor i<V(g 2 ) in eq. (|3.2| ) can be 
derived by studying the function U(p 2 ,p' 2 , q 2 ) in eq. (|3.4j ) and using in eqs. ( |3.9| ), ( |3.11| ) the 
following expressions: 



H = -2f K f } 



M 2 B 



m b M B + M K 



(Al) 



p(s,s,q ) 



2\l 



2A's - Au 



A 



6A'V + 2Ass' - QAA'su + A 2 u 2 



(A2) 



with A = s — ml, A' — s' — m 2 , u = s + s' — q 2 , A = u 2 — 4ss'. The coefficients of D = 3 
and D = 5 vacuum matrix elements are given by: 



4 = 



(A3) 



with r = p 2 — ml, r' = p' 2 — m 2 s . 

Also Fi(q 2 ) and -F (g 2 ) can be derived by equations analogous to 
relevant quantities for F\(q 2 ) are given by: 



H = Mb 



Ml 



(A4) 
1, O- The 

(A5) 



p(s,s', g 2 ^ 



m 6 [2A(tt - s') + A'(it - 4s)] + m s (Au - 2A's) 



+ 



2m b 



A 2 (3s , w-2ss / -w 2 ) + A /2 (3sw-6s 2 ) + 2AA'(3sw-2ss , - u 2 )] \ (A6) 



2vr' 



+ 



m; 



+ 



1 2m 2 — m s m b — 2q 2 



+ 



For Fo(q 2 ) the formulae read as follows: 



H = Mb^(M 2 b - M 2 K ) 



(A7) 
(A8) 

(A9) 



2A's - Am 



p(s, s', g 2 ) = ^= I A(m fe - m s ) + [m 6 (2(A - A') + 2s - u) + m s (w - 2s)] 



8^ 1 



A 



(A10) 
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m b (m b - m s ) 



rr 



(All) 



d$ = (m b - m s ) 



m b m 2 s m b {m b — m b m s + m 2 — q 2 ) 



+ 



m b - m s 
Qrr' 3 



+ 



mi 



+ 



2m h 



2r 3 r' 3r 2 r 



(A12) 



In the formulae for the coefficients of the non perturbative contributions, reported in this 
Appendix and in the following one, we have omitted all terms that vanish after the double 
Borel transform. 
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APPENDIX B: B 

The quantities appearing in the sum rule for the form factor Ti(q 2 ) in eq. (|4.2|) read as 
follows: 

M 2 

H = 2f K »M K *f B -* (Bl) 
m 6 



3 1 

p(s, s', q 2 ) = —/={^ - T-[m b m s (2A's + 2As' - u(A + A')) 



+ (s' - u)(2A's - Au) - s(2As' - A'u) - A' 2 s - A 2 s' + AA'm]} (B2) 



4 = -^=^ («>) 

m s 3m b + 2m s m 2 (m & + m s ) 
ml(m b + m s ) _ (m b + m s )[m b m s + 2(q 2 - mj - m 2 s )} 



For the form factor Ta(g ): 

M 2 

= —f'K*M K *f B — — (M| - M 2 -,) (B5) 
m 6 

p(s, s', q 2 ) = ^= |m fc m s (A' - A) + (As' - A's) + (u - 2s) + ^ ~ AA ' u) | (B6) 

(ra 6 - m,) [(m 6 + m s ) 2 - q 2 ] 
d * = 4^ ( B? ) 

3m b + m s (m fe - m s )m 2 [(mb + m s ) 2 - g 2 ] (m 6 - m s )m 2 [(m b + m s ) 2 - q 2 ] 
(m b - m s ) [(m b + m s ) 2 - g 2 ] [m b m s - 2{m\ + m 2 - g 2 )] 

+ ^^v 2 

q 2 (2m b — m s ) — 2{m b + m s )(m b — m 2 ) q 2 (5m b — 4m s ) — 4(m& + m s )(m 2 — m 2 ) 
24:7*7^'^ 24t >2 t 1 ' 

Form factor ^(g 2 ): 

M 2 2 

H = fx* Mx* f b — ~ (B9) 

JK KJB m b M B + M K * y 1 
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P = —^[m b (2As' - A'u) + m s {2Ms - Au)} 



4 = -A 

rr 



1 m\ m 2 2(m 2 b + m 2 s - q 2 ) - m b m s 

"5 — - n , + — ; + — + 



2^2^/ 27""^7"^ 2t"7"^^ ^i^*2^i/2 

Form factor Ai(g 2 ): 

M 2 

if = f K *M K *f B —?-(M B + M K *) 



3 

P = 



8^ 



(m b A' + m s A) + ^(A' 2 s + AV - AA'u) 
A 



d 3 = ~^r,\- m l + ni 2 -q 2 + 2m b m s ] 



1 3m 2 + 9m b m s + 4m 2 — 4g 2 4m 2 + 6m b m s + 6m 2 — 4g 2 
5 ~ ~6r? H 12rV ' 24r 2 r' 2 

m 2 [(m b + m s ) 2 - q 2 ] + m 2 [(m 6 + m s ) 2 - g 2 ] 



—4m 4 ~ 6mfm s — 4m 2 m 2 — Qm b m? s — 4m A s — 4m 2 g 2 + 6m b m s q 2 + 8m 2 g 2 — 4g 4 
Form factor A 2 (g 2 ): 

-H" = fK*M K *f B - 



m b M B + M/f. 



3 

p = -^72 { m b(2As' - A'u) + m s (2A's - Am) 



2m b 



A' 2 (2ss' + u 2 - 3su) + 3A 2 (2s' 2 - s'u) + 2AA'(-3s'w + 2ss' + u 2 )]} 



* = 1 



2rr' 

1 m 2 m 2 2m 2 + 2m 2 — 2g 2 — m b m s 



^^/y^/yf ^^/y/yf'^i ^jy^t/y? ^2^2^/2 
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TABLES 

TABLE I. Wilson coefficients d(n) for A^L = 225 MeV, n = 5 GeV and m t = 174 GeL. 



NDR 


Ci 


-0.243 


c 2 


1.105 


c 3 


1.083 x 10~ 2 


c 4 


-2.514 x 10~ 2 


c 5 


7.266 x 10~ 3 


c 6 


-3.063 x lO" 2 


cv 


-0.312 


c 9 


4.193 


Cio 


-4.578 



HV 



TABLE II. Parameters of the form factors. The functional q 2 dependence is either polar: 
F(q 2 ) = F(0)/(1 - q 2 /Ml) or linear: F{q 2 ) = F(0)(1 + (3q 2 ). For the form factor F T see text. 



F(0) M P (GeV) [3 [GeV- 2 ) 

J\ 0.25 ±0.03 5 

F 0.25 ± 0.03 7 

F T -0.14 ±0.03 

V 0.47 ±0.03 5 

A x 0.37 ±0.03 -0.023 

A 2 0.40 ± 0.03 0.034 

A 0.30 ± 0.03 4.8 

Ti 0.19 ±0.03 5.3 

T 2 0.19 ±0.03 -0.02 

T 3 -0.7 0.005 
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FIGURE CAPTIONS 



Fig. 1 

Form factors Fi(q 2 ), F (q 2 ) and i<V(g 2 ) of the transition B — > K£ + £~. The curves refer 
to different sets of parameters: s = 33 GeV^ 2 and s' Q = 1.4 GeV" 2 (continuous line), s = 
33 Ge^ 2 and s' = 1.6 Gel/ 2 (dashed line), s = 36 Ge\/ 2 and s = 1.4 GeV 2 (dotted line), 
so = 36 GeV 2 and s = 1.6 GeV 2 (dashed-dotted line). The Borel parameters are fixed to 
M 2 = 8 GeV 2 , M' 2 = 2 GeV 2 . 

Fig. 2 

Form factors Ti(q 2 ), T2(q 2 ) and T^(q 2 ) of the transition Z? — > K*£ + £~. The curves refer 
to different sets of parameters: so = 33 GeV^ 2 and s' = 1.6 GeV" 2 (continuous line), so = 
33 Gel/ 2 and s' = 1.8 Ge^ 2 (dashed line), s = 36 Gel/ 2 and s = 1.6 GeV 2 (dotted line), 
s = 36 GeV 2 and s = 1.8 GeV 2 (dashed-dotted line). The Borel parameters are fixed to 
M 2 = 8 GeV 2 , M' 2 = 2 GeV 2 . 

Fig. 3 

Form factors V(q 2 ), Ai(q 2 ), A 2 (q 2 ) and A (q 2 ) of -B — * K*£ + £~ . The curves refer to the 
same set of parameters as in fig. 2. 

Fig. 4 

Momentum dependence of the ratio between rare and semileptonic form factors TZ = 
Fi(q 2 ) I 'F? w \q 2 ); F t IW are obtained from eqs.CTPl). 

Fig. 5 

Invariant mass squared distribution of the lepton pair for the decay B — > K£ + £~: the dashed 
line refers to the short distance contribution only. 

Fig. 6 

Forward-backward asymmetry in the decay B — > K*£ + £~; the dashed line refers to the short 
distance contribution only. 

Fig. 7 

Invariant mass squared distribution of the lepton pair for the decay B — ► K*£ + £~: the 
dashed line refers to the short distance contribution only. 
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